------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\swhitt\Desktop\Gender and Accountability Replication Files\SCT gender and accountabili
> ty replication data wide format log file.log
  log type:  text
 opened on:  30 Nov 2023, 15:52:49

. do "C:\Users\swhitt\Desktop\Gender and Accountability Replication Files\SCT gender and accountability repl
> ication data wide format do file.do"

. *Gender, Agency, and Accountability for ISIS Violence: Public Perspectives from Mosul, Iraq
. *Replication Instructions
. 
. *Vera Mironova and Sam Whitt
. 
. *Below are instructions for replicating all manuscript and online appendix tables and figures in STATA usi
> ng the wide format of the 
. *dataset “SCT gender and accountability replication data wide format.dta”. 
. 
. *Please contact Sam Whitt (swhitt@highpoint.edu) for questions regarding data replication. 
. 
. *Note: You may need to install STATA packages for cibar, catcibar  and iebaltab commands. 
. *Use findit with command name to identify and download the appropriate packets to install. To install comm
> ands enter the following:
. 
. *ssc install cibar
. *net install catcibar, from("https://aarondwolf.github.io/catcibar")
. *findit iebaltab
. 
. *Note: In addition, some graphs require additional formatting using filename.grec files with the graph pla
> y command. 
. *To format a graph, simply run the command to generate the graph in the do file in STATA, then open the “G
> raph Editor” in STATA and click on 
. *the GREEN “Play Recording” button, then select “Browse” to select the grec file from among Replication fi
> les. 
. *The name of the grec file is indicated in the note below the graph command in the do file for the specifi
> c graph you wish to format. 
. *This should automatically format the graph, which you may then save to a location of your choosing.
. 
. *Manuscript Text
. 
. *“A total of 401 subjects were sampled across 11 neighborhoods in Mosul (See Appendix Map for Mosul sampli
> ng locations). 
. *Among the IDPs, we purposefully sampled camps outside Mosul that were believed to house people who had su
> pported ISIS 
. *(3 out of 15 camps in Nineveh province).”
. 
. tab location if mosul==1

          Mosul |
  Neighborhood; |
       IDP camp |
       location |      Freq.     Percent        Cum.
----------------+-----------------------------------
       AL Wahda |         36        8.98        8.98
     Al Muthana |         30        7.48       16.46
  Al Nabi Sheet |         25        6.23       22.69
        Al Noor |         30        7.48       30.17
      Al Rafaai |         32        7.98       38.15
     Al Zanjili |         48       11.97       50.12
       Al Zohur |         33        8.23       58.35
         Dawasa |         40        9.98       68.33
      Mohandsin |         31        7.73       76.06
   Mosul Jadida |         55       13.72       89.78
         Summar |         41       10.22      100.00
----------------+-----------------------------------
          Total |        401      100.00

. tab location if mosul==0

          Mosul |
  Neighborhood; |
       IDP camp |
       location |      Freq.     Percent        Cum.
----------------+-----------------------------------
  As Salamyiah  |         52       26.67       26.67
    Hamam Alil  |         93       47.69       74.36
Qayyarah Jadah  |         50       25.64      100.00
----------------+-----------------------------------
          Total |        195      100.00

. 
. 
. *“Fieldwork in the first wave in Mosul was conducted between December 23-30, 2018, and fieldwork in the ca
> mps was completed between 
. *January 10-24, 2019.”
. 
. tab date if mosul==1

    date of |
  interview |      Freq.     Percent        Cum.
------------+-----------------------------------
  23dec2018 |         65       16.21       16.21
  24dec2018 |         80       19.95       36.16
  25dec2018 |         55       13.72       49.88
  26dec2018 |         61       15.21       65.09
  27dec2018 |         63       15.71       80.80
  29dec2018 |         36        8.98       89.78
  30dec2018 |         41       10.22      100.00
------------+-----------------------------------
      Total |        401      100.00

. tab date if mosul==0

    date of |
  interview |      Freq.     Percent        Cum.
------------+-----------------------------------
  10jan2019 |         42       21.54       21.54
  12jan2019 |         51       26.15       47.69
  17jan2019 |         50       25.64       73.33
  24jan2019 |         52       26.67      100.00
------------+-----------------------------------
      Total |        195      100.00

. 
. *“In general, Mosul residents are more educated, have higher incomes, and tend to come from more professio
> nal backgrounds of employment or 
. *are students. Camp members are conversely of lower education, income, and tend to have manual labor backg
> rounds, which anecdotally fits 
. *into the public perception in Mosul of the typical ISIS demographic.”
. 
. iebaltab female age education income professional laborer unemployed student moved, groupvar(mosul) savexl
> s(mosulbalance)

    Balance table saved in Excel format to: mosulbalance.xlsx

. 
. *“Despite these limitations, 1 in 5 camp members acknowledged that they worked directly for ISIS or fought
>  for ISIS. 
. *Over half indicate that a family member worked for ISIS, while 1 in 5 also have a family member who fough
> t for ISIS, 
. *and nearly two-thirds have family members who have been put on trial for their role in ISIS. 
. *In total, 91% of camp IDPs indicated some form of ISIS connection or affiliation.”
. 
. sum workedisis-isisties

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
  workedisis |        195    .2102564    .4085396          0          1
  spouseisis |        195    .0358974    .1865133          0          1
     famisis |        195    .5487179    .4989018          0          1
  foughtisis |        195    .1692308    .3759208          0          1
spousefoug~s |        195    .1128205    .3171876          0          1
-------------+---------------------------------------------------------
famfoughti~s |        195          .2    .4010296          0          1
famtrialisis |        195    .6410256    .4809344          0          1
    isisties |        195    .9128205     .282824          0          1

. 
. *“Consistent with our sampling strategy, we find that nearly two-thirds of camp members believe that “almo
> st all” or 
. *“most people in Mosul supported ISIS” (63%), while very few respondents we sampled in Mosul (15%) stated 
> this was the case.”
. 
. tab supportisis if mosul==0

            9. Support for ISIS |      Freq.     Percent        Cum.
--------------------------------+-----------------------------------
                Almost everyone |         12        6.15        6.15
                    Most people |        111       56.92       63.08
Some people, but not a majority |         71       36.41       99.49
           Almost no one at all |          1        0.51      100.00
--------------------------------+-----------------------------------
                          Total |        195      100.00

. tab supportisis if mosul==1

            9. Support for ISIS |      Freq.     Percent        Cum.
--------------------------------+-----------------------------------
                Almost everyone |          9        2.24        2.24
                    Most people |         52       12.97       15.21
Some people, but not a majority |        326       81.30       96.51
           Almost no one at all |         14        3.49      100.00
--------------------------------+-----------------------------------
                          Total |        401      100.00

. 
. *“To measure beliefs about women’s agency, we utilize a series of questions related to the role played by 
> ISIS men and women in the organization….
. *Factor analysis indicates that responses to these items are also highly inter-correlated (SI Table 2). 
. *We normalize and score these items into a common index called agency affirmation ranging from 0=strong de
> nial of women’s agency to 1=strong affirmation.”
. 
. *Agency Affirmation Variable
. 
. factor revwomenradical revwomendanger revwomenorg revwomenhelpfight womenfollow menlead
(obs=562)

Factor analysis/correlation                      Number of obs    =        562
    Method: principal factors                    Retained factors =          3
    Rotation: (unrotated)                        Number of params =         15

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      2.88198      2.40524            0.9762       0.9762
        Factor2  |      0.47675      0.39805            0.1615       1.1376
        Factor3  |      0.07869      0.20913            0.0267       1.1643
        Factor4  |     -0.13043      0.02324           -0.0442       1.1201
        Factor5  |     -0.15367      0.04729           -0.0521       1.0681
        Factor6  |     -0.20096            .           -0.0681       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(15) = 1436.26 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    -----------------------------------------------------------
        Variable |  Factor1   Factor2   Factor3 |   Uniqueness 
    -------------+------------------------------+--------------
    revwomenra~l |   0.7424    0.0397   -0.1602 |      0.4217  
    revwomenda~r |   0.7725    0.1370   -0.1294 |      0.3678  
     revwomenorg |   0.7196    0.2862    0.0766 |      0.3944  
    revwomenhe~t |   0.5463    0.3045    0.1520 |      0.5858  
     womenfollow |   0.6628   -0.3908    0.0473 |      0.4057  
         menlead |   0.6918   -0.3592    0.0714 |      0.3872  
    -----------------------------------------------------------

. 
. *predict agencyaffirm
. *norm agencyaffirm, method(mmx)
. 
. *“To measure beliefs about the judicial system’s affirmation or denial of women’s agency, we ask responden
> ts whether they agree or disagree 
. *that different groups are being sufficiently punished for their role in ISIS on a four-point scale rangin
> g from strongly agree to strongly 
. *disagree…We then create a normalized index based on responses to questions specifically about women calle
> d judicial denial of agency, 
. *ranging from 0=strong judicial affirmation (agreement that women are being sufficiently punished) to 
. *1=strong judicial denial of agency (disagreement that women are being sufficiently punished).”
. 
. *Denial of Judicial Agency Variable
. 
. factor leaderswives fighterwives femaleworkers
(obs=591)

Factor analysis/correlation                      Number of obs    =        591
    Method: principal factors                    Retained factors =          1
    Rotation: (unrotated)                        Number of params =          3

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      1.69783      1.75531            1.1526       1.1526
        Factor2  |     -0.05749      0.10986           -0.0390       1.1136
        Factor3  |     -0.16734            .           -0.1136       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(3)  =  699.76 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    ---------------------------------------
        Variable |  Factor1 |   Uniqueness 
    -------------+----------+--------------
    leaderswives |   0.8162 |      0.3338  
    fighterwives |   0.8388 |      0.2963  
    femalewor~rs |   0.5727 |      0.6720  
    ---------------------------------------

. 
. *predict judicialdenial
. *norm judicialdenial, method(mmx)
. 
. *“There are two observations for each individual in the study: a punishment of a male ISIS worker and a fe
> male ISIS worker, 
. *and standard errors are clustered by respondent.”
. 
. *To generate the dependent variable (see long format), stack and combine responses to variables “nasimopp”
>  and “saraworker” for a 
. *within-subject comparison for male and female isis workers.
. 
. *“Among all respondents, punishment preferences for women workers averaged 2.77 (sd.=1.43) compared to 2.3
> 4 for men (sd=1.21), a significant difference based on paired t-tests (t=6.42, p<0.0000).”
. 
. ttest maleworker = femaleworker

Paired t test
------------------------------------------------------------------------------
Variable |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
malewo~r |     344    2.351744    .0652552    1.210304    2.223393    2.480095
female~r |     344    2.773256    .0771506    1.430931    2.621508    2.925004
---------+--------------------------------------------------------------------
    diff |     344   -.4215116    .0656077    1.216841   -.5505556   -.2924676
------------------------------------------------------------------------------
     mean(diff) = mean(maleworker - femaleworker)                 t =  -6.4247
 Ho: mean(diff) = 0                              degrees of freedom =      343

 Ha: mean(diff) < 0           Ha: mean(diff) != 0           Ha: mean(diff) > 0
 Pr(T < t) = 0.0000         Pr(|T| > |t|) = 0.0000          Pr(T > t) = 1.0000

. 
. *“First, victimization by ISIS does not predict more severe punishments.” 
. 
. *ISIS victimization index
. 
. factor punishedis- womenabusedis
(obs=596)

Factor analysis/correlation                      Number of obs    =        596
    Method: principal factors                    Retained factors =          5
    Rotation: (unrotated)                        Number of params =         35

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      1.29834      0.54391            0.7524       0.7524
        Factor2  |      0.75444      0.35373            0.4372       1.1896
        Factor3  |      0.40071      0.35703            0.2322       1.4219
        Factor4  |      0.04368      0.03494            0.0253       1.4472
        Factor5  |      0.00874      0.10054            0.0051       1.4522
        Factor6  |     -0.09180      0.07987           -0.0532       1.3990
        Factor7  |     -0.17166      0.07695           -0.0995       1.2996
        Factor8  |     -0.24861      0.01967           -0.1441       1.1555
        Factor9  |     -0.26828            .           -0.1555       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(36) =  633.33 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    -------------------------------------------------------------------------------
        Variable |  Factor1   Factor2   Factor3   Factor4   Factor5 |   Uniqueness 
    -------------+--------------------------------------------------+--------------
      punishedis |   0.3398    0.4775   -0.1422   -0.0247    0.0141 |      0.6355  
    familypuni~s |   0.2723    0.4101   -0.1536    0.0832   -0.0280 |      0.7264  
       injuredis |  -0.1843   -0.1268    0.2567    0.1292    0.0025 |      0.8674  
    familyinju~s |   0.2058    0.0537    0.3332    0.0159   -0.0465 |      0.8413  
    familykill~s |   0.2273    0.0640    0.3169   -0.0575    0.0500 |      0.8380  
    imprisonedis |   0.3036    0.1660    0.0723    0.1040    0.0377 |      0.8628  
      fleehomeis |   0.6125   -0.3829   -0.0616   -0.0037   -0.0267 |      0.4738  
     hometakenis |   0.6455   -0.3271   -0.0992    0.0177    0.0209 |      0.4658  
    womenabuse~s |   0.3111    0.2325    0.2463   -0.0689   -0.0227 |      0.7833  
    -------------------------------------------------------------------------------

. *predict isisvictim
. 
. *Manuscript Tables and Figures
. 
. 
. *Figure 1. ISIS punishment preferences (no gender prime)
. *see long format
. 
. 
. *Note additional formatting requires the "Figure 1 long formatting.grec" file with the command graph play 
> "Figure 1 long formatting.grec"
. 
. *Figure 2. ISIS punishment preferences (gender primes)
. *see long format
. 
. *Figure 3. Women's agency and accountability
. 
. *Wide format
. 
. cibar agencyaffirm, over1(mosul) 

. graph save g1, replace
(file g1.gph saved)

. cibar judicialdenial, over1(mosul)

. graph save g2, replace
(file g2.gph saved)

. graph combine "g1.gph" "g2.gph"

. 
. *Note additional formatting requires the "Figure 3 wide formatting.grec" file with the command graph play 
> "Figure 3 wide formatting.grec"
. 
. *Table 2. Women’s Agency for Violence and Accountability (Ordered Probit Regression)
. *see long format
. 
. *Appendix Replication
. 
. *Mosul Sampling Locations
. *Wide format
. 
. tab location if mosul==1

          Mosul |
  Neighborhood; |
       IDP camp |
       location |      Freq.     Percent        Cum.
----------------+-----------------------------------
       AL Wahda |         36        8.98        8.98
     Al Muthana |         30        7.48       16.46
  Al Nabi Sheet |         25        6.23       22.69
        Al Noor |         30        7.48       30.17
      Al Rafaai |         32        7.98       38.15
     Al Zanjili |         48       11.97       50.12
       Al Zohur |         33        8.23       58.35
         Dawasa |         40        9.98       68.33
      Mohandsin |         31        7.73       76.06
   Mosul Jadida |         55       13.72       89.78
         Summar |         41       10.22      100.00
----------------+-----------------------------------
          Total |        401      100.00

. tab location if mosul==0

          Mosul |
  Neighborhood; |
       IDP camp |
       location |      Freq.     Percent        Cum.
----------------+-----------------------------------
  As Salamyiah  |         52       26.67       26.67
    Hamam Alil  |         93       47.69       74.36
Qayyarah Jadah  |         50       25.64      100.00
----------------+-----------------------------------
          Total |        195      100.00

. 
. *Summary of Variables in Manuscript Analysis
. *Wide format
. 
. sum mosul maleworker femaleworker agencyaffirm judicialdenial isisvictim female age education income profe
> ssional laborer student unemployed moved

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
       mosul |        596    .6728188    .4695781          0          1
  maleworker |        344    2.351744    1.210304          1          5
femaleworker |        345    2.771014    1.429456          1          5
agencyaffirm |        562    .5963696    .2510677          0          1
judiciald~al |        591    .5284018    .2967733          0          1
-------------+---------------------------------------------------------
  isisvictim |        596    4.11e-09    .8120216   -.833234   3.157621
      female |        596     .216443    .4121655          0          1
         age |        596    34.71309    12.75456         18         72
   education |        595    2.606723    .9137678          1          4
      income |        596    2.130872    .9310624          1          4
-------------+---------------------------------------------------------
professional |        596    .2567114    .4371857          0          1
     laborer |        596    .3305369    .4708016          0          1
     student |        596     .238255    .4263738          0          1
  unemployed |        596    .1694631    .3754756          0          1
       moved |        596    .1795302     .384118          0          1

. 
. *Figure 1 Box-Whisker Plots of Punishing ISIS from Manuscript Figure 1
. *Wide format
. 
. graph box dleaders-dtaxpayers, by(mosul) nolabel

. 
. *Note additional formatting requires the "SA Figure 1a wide formatting.grec" file with the command graph p
> lay "SA Figure 1a wide formatting.grec"
. 
. catcibar dleaders- dtaxpayers, by(mosul)

. *Note additional formatting requires the "SA Figure 1b wide formatting.grec" file with the command graph p
> lay "SA Figure 1b wide formatting.grec"
. 
. *(Alternative Version) see Long-format
. 
. *Figure 2 Box-Whisker Plots of Vignette Punishments from Manuscript Figure 2.
. *see also long format
. *Wide format
. 
. catcibar loyalfighter loyalworker wifeworker wifelowhigh , over(treatment) by(mosul)

. *Note additional formatting requires the "SA Figure 2b wide formatting.grec" file with the command graph p
> lay "SA Figure 2b wide formatting.grec"
. 
. *Figure 3. Perceptions of Judicial Affirmation/Denial of Women’s Agency
. *see long format
. 
. *Table 1. Factor Analysis of Judicial Affirmation/Denial Index
. *Wide format
. 
. factor punishedleaders- femaleworkers
(obs=563)

Factor analysis/correlation                      Number of obs    =        563
    Method: principal factors                    Retained factors =          3
    Rotation: (unrotated)                        Number of params =         18

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      2.58423      1.85525            0.9005       0.9005
        Factor2  |      0.72899      0.62819            0.2540       1.1545
        Factor3  |      0.10079      0.14352            0.0351       1.1896
        Factor4  |     -0.04272      0.06425           -0.0149       1.1748
        Factor5  |     -0.10697      0.05017           -0.0373       1.1375
        Factor6  |     -0.15714      0.08027           -0.0548       1.0827
        Factor7  |     -0.23740            .           -0.0827       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(21) = 1324.13 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    -----------------------------------------------------------
        Variable |  Factor1   Factor2   Factor3 |   Uniqueness 
    -------------+------------------------------+--------------
    punishedle~s |   0.4088    0.4955    0.1245 |      0.5719  
    punishedfi~s |   0.4317    0.2389   -0.1670 |      0.7287  
     punishedciv |   0.6460    0.0240   -0.1252 |      0.5664  
      punishedff |   0.5085    0.4628    0.0189 |      0.5269  
    leaderswives |   0.8058   -0.2437   -0.0509 |      0.2887  
    fighterwives |   0.7644   -0.3605    0.0303 |      0.2848  
    femalewor~rs |   0.5668   -0.1492    0.1946 |      0.6186  
    -----------------------------------------------------------

. factor leaderswives fighterwives femaleworkers
(obs=591)

Factor analysis/correlation                      Number of obs    =        591
    Method: principal factors                    Retained factors =          1
    Rotation: (unrotated)                        Number of params =          3

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      1.69783      1.75531            1.1526       1.1526
        Factor2  |     -0.05749      0.10986           -0.0390       1.1136
        Factor3  |     -0.16734            .           -0.1136       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(3)  =  699.76 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    ---------------------------------------
        Variable |  Factor1 |   Uniqueness 
    -------------+----------+--------------
    leaderswives |   0.8162 |      0.3338  
    fighterwives |   0.8388 |      0.2963  
    femalewor~rs |   0.5727 |      0.6720  
    ---------------------------------------

. 
. *Figure 4. Perceptions of Women’s Agency in ISIS
. *see long format
. 
. *Table 2. Factor Analysis of Agency Affirmation Index
. *Wide format
. 
. factor revwomenradical revwomendanger revwomenorg revwomenhelpfight womenfollow menlead
(obs=562)

Factor analysis/correlation                      Number of obs    =        562
    Method: principal factors                    Retained factors =          3
    Rotation: (unrotated)                        Number of params =         15

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      2.88198      2.40524            0.9762       0.9762
        Factor2  |      0.47675      0.39805            0.1615       1.1376
        Factor3  |      0.07869      0.20913            0.0267       1.1643
        Factor4  |     -0.13043      0.02324           -0.0442       1.1201
        Factor5  |     -0.15367      0.04729           -0.0521       1.0681
        Factor6  |     -0.20096            .           -0.0681       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(15) = 1436.26 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    -----------------------------------------------------------
        Variable |  Factor1   Factor2   Factor3 |   Uniqueness 
    -------------+------------------------------+--------------
    revwomenra~l |   0.7424    0.0397   -0.1602 |      0.4217  
    revwomenda~r |   0.7725    0.1370   -0.1294 |      0.3678  
     revwomenorg |   0.7196    0.2862    0.0766 |      0.3944  
    revwomenhe~t |   0.5463    0.3045    0.1520 |      0.5858  
     womenfollow |   0.6628   -0.3908    0.0473 |      0.4057  
         menlead |   0.6918   -0.3592    0.0714 |      0.3872  
    -----------------------------------------------------------

. 
. *Figure 5. Perceptions of ISIS Victimization by Location
. *Wide format
. 
. catcibar punishedis- womenabusedis, over(mosul)

. 
. *Table 3. Factor Analysis of Victimization Index
. *Wide format
. 
. factor punishedis- womenabusedis
(obs=596)

Factor analysis/correlation                      Number of obs    =        596
    Method: principal factors                    Retained factors =          5
    Rotation: (unrotated)                        Number of params =         35

    --------------------------------------------------------------------------
         Factor  |   Eigenvalue   Difference        Proportion   Cumulative
    -------------+------------------------------------------------------------
        Factor1  |      1.29834      0.54391            0.7524       0.7524
        Factor2  |      0.75444      0.35373            0.4372       1.1896
        Factor3  |      0.40071      0.35703            0.2322       1.4219
        Factor4  |      0.04368      0.03494            0.0253       1.4472
        Factor5  |      0.00874      0.10054            0.0051       1.4522
        Factor6  |     -0.09180      0.07987           -0.0532       1.3990
        Factor7  |     -0.17166      0.07695           -0.0995       1.2996
        Factor8  |     -0.24861      0.01967           -0.1441       1.1555
        Factor9  |     -0.26828            .           -0.1555       1.0000
    --------------------------------------------------------------------------
    LR test: independent vs. saturated:  chi2(36) =  633.33 Prob>chi2 = 0.0000

Factor loadings (pattern matrix) and unique variances

    -------------------------------------------------------------------------------
        Variable |  Factor1   Factor2   Factor3   Factor4   Factor5 |   Uniqueness 
    -------------+--------------------------------------------------+--------------
      punishedis |   0.3398    0.4775   -0.1422   -0.0247    0.0141 |      0.6355  
    familypuni~s |   0.2723    0.4101   -0.1536    0.0832   -0.0280 |      0.7264  
       injuredis |  -0.1843   -0.1268    0.2567    0.1292    0.0025 |      0.8674  
    familyinju~s |   0.2058    0.0537    0.3332    0.0159   -0.0465 |      0.8413  
    familykill~s |   0.2273    0.0640    0.3169   -0.0575    0.0500 |      0.8380  
    imprisonedis |   0.3036    0.1660    0.0723    0.1040    0.0377 |      0.8628  
      fleehomeis |   0.6125   -0.3829   -0.0616   -0.0037   -0.0267 |      0.4738  
     hometakenis |   0.6455   -0.3271   -0.0992    0.0177    0.0209 |      0.4658  
    womenabuse~s |   0.3111    0.2325    0.2463   -0.0689   -0.0227 |      0.7833  
    -------------------------------------------------------------------------------

. 
. *Further Analysis of Gender and ISIS Punishment Preferences
. *see long format
. 
. *Table 4. Punishment of ISIS (OLS Regression)
. *see long format
. 
. *Table 5. Punishment of ISIS (Ordered Probit Regression)
. *see long format
. 
. *Table 6. Punishment of ISIS by Treatment Groups (OLS regression)
. *Wide format
. 
. reg loyalworker treatment mosul judicialdenialall isisvictim female age education income professional labo
> rer unemployed moved, robust

Linear regression                               Number of obs     =        561
                                                F(12, 548)        =      19.56
                                                Prob > F          =     0.0000
                                                R-squared         =     0.2365
                                                Root MSE          =      1.101

-----------------------------------------------------------------------------------
                  |               Robust
      loyalworker |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
        treatment |  -.1202533   .1013109    -1.19   0.236    -.3192585    .0787519
            mosul |   .8466913   .1432771     5.91   0.000     .5652518    1.128131
judicialdenialall |   .4708819   .2360265     2.00   0.047     .0072544    .9345094
       isisvictim |   .0466411   .0628059     0.74   0.458    -.0767287    .1700109
           female |   -.310697   .1108814    -2.80   0.005    -.5285015   -.0928925
              age |   .0016677   .0046271     0.36   0.719    -.0074214    .0107568
        education |  -.0148248   .0725247    -0.20   0.838    -.1572853    .1276356
           income |   .1815904   .0726254     2.50   0.013      .038932    .3242487
     professional |   .0960977   .1628636     0.59   0.555    -.2238156     .416011
          laborer |  -.0166849   .1590478    -0.10   0.916    -.3291028     .295733
       unemployed |   .1866572   .1747798     1.07   0.286    -.1566633    .5299776
            moved |   .0642072   .1235159     0.52   0.603    -.1784153    .3068298
            _cons |   1.521275   .3294071     4.62   0.000     .8742195     2.16833
-----------------------------------------------------------------------------------

. 
. reg loyalfighter treatment mosul judicialdenialall isisvictim female age education income professional lab
> orer unemployed moved, robust

Linear regression                               Number of obs     =        562
                                                F(12, 549)        =      91.99
                                                Prob > F          =     0.0000
                                                R-squared         =     0.6180
                                                Root MSE          =     .81236

-----------------------------------------------------------------------------------
                  |               Robust
     loyalfighter |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
------------------+----------------------------------------------------------------
        treatment |    .062145   .0724182     0.86   0.391    -.0801057    .2043958
            mosul |   2.094921   .1085081    19.31   0.000     1.881779    2.308063
judicialdenialall |   .2031473   .1497458     1.36   0.175    -.0909975    .4972922
       isisvictim |  -.1080814   .0522518    -2.07   0.039    -.2107193   -.0054435
           female |  -.2159175   .0971121    -2.22   0.027    -.4066743   -.0251606
              age |  -.0009425    .003142    -0.30   0.764    -.0071142    .0052293
        education |   .0554283   .0471383     1.18   0.240    -.0371652    .1480219
           income |    .087824    .050512     1.74   0.083    -.0113965    .1870445
     professional |   .2686208   .1221859     2.20   0.028     .0286118    .5086298
          laborer |   .2409145   .1170745     2.06   0.040     .0109457    .4708832
       unemployed |   .0954867   .1204001     0.79   0.428    -.1410144    .3319879
            moved |   .0191886   .0905365     0.21   0.832    -.1586518     .197029
            _cons |   1.756834   .2299943     7.64   0.000     1.305058    2.208611
-----------------------------------------------------------------------------------

. 
. reg wifeworker treatment mosul judicialdenial agencyaffirm isisvictim female age education income professi
> onal laborer unemployed moved, robust

Linear regression                               Number of obs     =        556
                                                F(13, 542)        =      29.56
                                                Prob > F          =     0.0000
                                                R-squared         =     0.3331
                                                Root MSE          =     1.1163

--------------------------------------------------------------------------------
               |               Robust
    wifeworker |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
     treatment |   .3981426   .1055458     3.77   0.000     .1908136    .6054716
         mosul |   .9135491   .1544028     5.92   0.000     .6102478     1.21685
judicialdenial |   .3134052   .1751113     1.79   0.074    -.0305748    .6573852
  agencyaffirm |   1.121978   .2653531     4.23   0.000     .6007317    1.643225
    isisvictim |   .0801109   .0677475     1.18   0.238    -.0529689    .2131906
        female |  -.0567379   .1171041    -0.48   0.628    -.2867714    .1732957
           age |   .0009249   .0050857     0.18   0.856    -.0090652     .010915
     education |  -.0538253   .0707847    -0.76   0.447    -.1928712    .0852206
        income |   .1275584   .0750756     1.70   0.090    -.0199163    .2750332
  professional |   .1628465   .1707562     0.95   0.341    -.1725785    .4982715
       laborer |  -.1692807    .160196    -1.06   0.291    -.4839617    .1454003
    unemployed |  -.0840884   .1546702    -0.54   0.587    -.3879148     .219738
         moved |   .1179501   .1177919     1.00   0.317    -.1134344    .3493346
         _cons |   .6132663   .3186336     1.92   0.055    -.0126417    1.239174
--------------------------------------------------------------------------------

. 
. reg wifelowhigh treatment mosul judicialdenial agencyaffirm isisvictim female age education income profess
> ional laborer unemployed moved, robust

Linear regression                               Number of obs     =        548
                                                F(13, 534)        =      29.57
                                                Prob > F          =     0.0000
                                                R-squared         =     0.3092
                                                Root MSE          =     1.1691

--------------------------------------------------------------------------------
               |               Robust
   wifelowhigh |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
---------------+----------------------------------------------------------------
     treatment |   .3304846   .1122053     2.95   0.003     .1100665    .5509026
         mosul |   .9912808   .1884355     5.26   0.000     .6211151    1.361447
judicialdenial |   .6089888   .1955629     3.11   0.002     .2248219    .9931558
  agencyaffirm |   .9777026   .3081297     3.17   0.002     .3724076    1.582998
    isisvictim |   .1445989   .0747585     1.93   0.054     -.002258    .2914558
        female |  -.1938089   .1270288    -1.53   0.128    -.4433463    .0557284
           age |   -.007807   .0054417    -1.43   0.152    -.0184968    .0028828
     education |   -.059333   .0736721    -0.81   0.421    -.2040556    .0853896
        income |   .0241089   .0715637     0.34   0.736    -.1164721    .1646899
  professional |   .2953673   .1751814     1.69   0.092    -.0487619    .6394965
       laborer |  -.0911475   .1696885    -0.54   0.591    -.4244863    .2421914
    unemployed |   .0028293   .1636198     0.02   0.986    -.3185882    .3242468
         moved |   .1539403   .1290841     1.19   0.234    -.0996345    .4075152
         _cons |   .9648346   .3347135     2.88   0.004      .307318    1.622351
--------------------------------------------------------------------------------

. 
. *Table  7. Summary of Additional Variables used in SI Tables 4-6
. *see long format
. *wide format
. 
. sum loyalworker loyalfighter wifeworker wifelowhigh

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
 loyalworker |        593    2.543002    1.236517          1          5
loyalfighter |        596    3.781879    1.325664          1          5
  wifeworker |        596    2.805369    1.356852          1          5
 wifelowhigh |        588    2.727891    1.395089          1          5

. 
. *Manuscript Table 2 Robustness Checks
. *see long format
. 
. *Manuscript Table 2 Power Calculations
. *Wide format
. 
. pwcorr maleworker femaleworker

             | malewo~r female~r
-------------+------------------
  maleworker |   1.0000 
femaleworker |   0.5866   1.0000 

. 
. power pairedmeans 2.35 2.77, corr(0.5866) sd1(1.21) sd2(1.43) power(0.80 0.90 0.95 0.99) graph

. 
. log close
      name:  <unnamed>
       log:  C:\Users\swhitt\Desktop\Gender and Accountability Replication Files\SCT gender and accountabili
> ty replication data wide format log file.log
  log type:  text
 closed on:  30 Nov 2023, 15:53:06
------------------------------------------------------------------------------------------------------------
